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Abstract 

We apply the chiral color dielectric model to low- and medium-energy scatter- 
ing within the coupled vrN and vrA system. Dynamic baryon states in which 
quarks are confined by the scalar color-dielectric field are constructed in Fock 
space. Spurious motion of the center of mass is eliminated by construct- 
ing momentum eigenstates via a Peierls-Yoccoz projection. The relativistic 
Lippmann-Schwinger equation is solved for the complex energies of the T ma- 
trix poles, and the pole positions of the nucleon and delta are used to fix the 
few parameters of the model. The S- and P-wave phase shifts, the bare N 
and A masses, the renormalized vrNN and ttNA coupling constants, and the 
vrNN and vrNA vertex functions are predicted. 
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I. INTRODUCTION 



From their first introduction, quark models have been extensively used to study baryon 
spectra. Some models have the quarks bound by a potential [Q] , others have them confined 
to a bag 0], and others have them bound as solitons 0. These models are fairly successful 
at fitting the masses and static properties of the baryon octet and decuplet and even some 
higher excited states. In most quark models, the baryon masses are identified with the 
bound state energies of quarks within some potential well. This means that the pion-baryon 
interactions do not directly affect the values of the baryon masses, and that the baryons are 
stable since decay channels are not coupled in. Phenomenologically, the masses of excited 
baryons are determined as the resonance energies in the phase-shift analysis of pion-nucleon 
scattering. These experimental masses naturally include mass shifts due to the pion-baryon 
interactions as well as widths due to the finite lifetimes of the states. Accordingly, there are 
uncertainties in a direct comparison of the experimental masses to the bound-state energies 
of quarks bound in a potential well. 

Chiral models, such as the cloudy bag models 0, contain pions in additions to quarks, 
and therefore also contain elementary pion-quark interactions. The calculation of masses 
with these models consequently include shifts due to the pion-baryon interactions. However, 
these models usually use perturbation theory to calculate the widths of the excited baryons, 
and this means that higher-order terms are ignored and that unitarity in the pion-nucleon 
channel is not ensured. Ideally, one should determine the masses and widths of excited 
baryons by calculating the pion-baryon scattering matrix for a given model, and then finding 
the complex energies at which the scattering matrix has poles. This approach was used by 
Thomas, Theberge and Miller |^ for the N-A system, and later extended to other partial 
waves and to other mesons [0. 

Bag models suffer from several limitations. First, the use of an artificial bag to confine 
quarks leads to wave functions which have sharp discontinuities at the bag boundary. This 
leads to pion-baryon form factors (and ultimately potentials) which have unphysical oscil- 
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lations and, presumably, limited ability to describe the scattering data. In addition, bag 
calculations treat the bare baryon masses and the bag radius as independent parameters 
even though we know this is not correct since for massless quarks the baryon masses are 
inversely proportional to the bag radius 0. 

In the present work the chiral color dielectric model (CDM) is employed to overcome 
some of the limitations of bag models. The CDM was first introduced phenomenologically to 
simulate the absolute confinement of quantum chromodynamics (QCD). A connection with 
the more fundamental theory was established by Nielson and Patkos [§] who found that 
an effective color dielectric field does indeed arise naturally in lattice QCD after defining 
"coarse grained" effective field variables. The color dielectric field takes into account the long 
distance behavior of the QCD vacuum and produces a natural confinement of the quarks 
within baryons. More recently, Krein et al. have solved the Schwinger-Dyson equation 
and shown how confinement and a pseudoscalar meson (Goldstone boson) arise if chiral 
symmetry is dynamically broken. 

In some recent CDM models, chiral symmetry is restored via a nonlinear realization 
which introduces an elementary pion field as a Goldstone boson. Even though the mesons 
in these models are elementary, the baryon structure and the resulting relations of coupling 



constants are determined by the quark model |[10|JTl|| . In application, the CDM is similar to 



the cloudy bag model yet somewhat more microscopic in its dynamic mechanism for 

quark confinement and its treatment of baryon recoil within a nontopological-soliton solution 
to the field equations. Accordingly, when in the present paper we compare models which 
differ mainly in their treatment of quark binding, we hope to determine how important are 
improved (or maybe just different) quark dynamics when calculating two-body scattering, 
and how might the differences in the models manifest themselves. 

Our calculation proceeds as follows. First we solve the quark and dielectric field equations 
of CDM in the mean field approximation (MFA). We then quantize the classical dielectric 
field using the coherent state, and improve upon the MFA by including one-gluon exchange 
corrections perturbatively. We include the recoil of the baryon via the prescription of Peierls- 



Yoccoz projection [|T2|. This, presumably, is a more important correction to scattering 
where large momentum transfers can occur than for static properties where the momentum 
transfers are small. Once we have bare baryons constructed from quarks within a binding 
field, we expand the quark-pion interaction in powers of the pion field and calculate the 
transition matrix elements for baryon coupling with multiple pions. From these we identify 
the pion-nucleon potential, and then dress the vertices by substituting the potentials into 
the relativistic Lippmann-Schwinger equation 0: 

Here the superscripts LJI indicate the orbital angular momenta, the total angular momenta, 
and the total isospin, and the subscripts /3 and a indicate the vrN or vrA channels. We 
adjust the parameters of color dielectric model so that the poles of T matrix elements occur 
at the nucleon mass and at the A's complex mass. We then use the model to predict low- 
and medium-energy vrN scattering in 5* and P waves and the ttNN and the ttNA coupling 
constants. 

II. CHIRAL COLOR DIELECTRIC MODEL 



The CDM Lagrangian we choose to use is: |10| 

C = q 



X 2 



cr. 
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-^ 2 2 ' 

<X) = X\ (3) 

D^4> = (f)d^(l) + f sin jd^(f), (4) 
F;, = d,At - d^Al - gfabcAlAl. (5) 
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Here q, (j), and x are the quark, pion, gluon, and dielectric fields respectively, k{x) is the 
dielectric coefficient (other CDM models may have different functional forms for ^Ilu 
is the gluon field tensor, r is the isospin operator of the pion, Xa and fabc are Gell-Mann 
matrices and SU(3) structrue constants, rrig and are the quark and pion masses, and D^(f) 
is the covariant derivative of the pion field. Arrows over quantities denote their isovector 
nature. The form @ of the Lagrangian becomes invariant under chiral transformations 
when the pion mass is set equal to zero. 

An important, yet phenomenological, part of the CDM Lagrangian is the dielectric 
self-interaction field U{x)- We express it in terms of the bag constant B, a shape parameter 
a, and the dielectric field x- 



U{x) = Bax' 



(6) 



This self-interaction field has an absolute minimum when the dielectric field x = and 
a > 6, and has a local minimum when % = 1. Of particular importance is the behavior of 
the Lagrangian (Q) in the limit of vanishing dielectric field, x 0. In this limit, the gluon 
kinetic energy term vanishes and the quark effective mass mq/x becomes infinite. 

Accordingly, we have the distinguishing characteristic of the CDM: the quark and gluon 
fields are confined within the baryon to the region of non- vanishing dielectric field x- This 
contrasts with the bag model in which confinement is externally imposed by placing the 
quarks within an infinite square well. 

As we see from the form of the Lagrangian @ with the dielectric self-interaction field 
(^), the CDM has five parameters: the quark mass nig, the strong coupling constant g (or 
equivalently = g^/Air), the pion decay constant /, the glueball (or dielectric field) mass 



irigb = \^2Ba/av, and the bag constant B. Although we do adjust some of these parameters 
to experimental data, their values are far from arbitrary. For example, the glueball mass 
nigi, essentially sets the scale of energy, and the values of the coupling and decay constants 
must be close to the accepted values for the model to be considered realistic. 
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A. Pseudo-vector Coupling 



In order for the Weinberg- Tomozawa relation to appear explicitly in S-wave vrN scattering 
at the tree level (and thus guarantee that the Born approximation scattering lengths are 
reasonable), we perform a unitary transformation ^ on the quark fields, q —>■ exp[i75r • 
(f)/{2f)]q. This transforms the original CDM Lagrangian (|^) to 



C = q 



TTl 1 

X ^ J 



(7) 



where the covariant derivative on the quark fields are defined as 



= daq 



cos 



l]r-d)xdu 



(8) 



These new quark fields are dressed by the pions. We next expand the Lagrangian in powers 
of the inverse of /, and retain terms up to order 1//^. We separate off terms linear and 
quadratic in 1// and identify them as the interaction Hamiltonian: 



(9) 



where 



H2n = ^ j d^xq-ff'Tq ■ (j) x d^, 



(10) 
(11) 



These elementary vertices are illustrated in Figure [l|. You will notice that the one-pion 
vertex Hit^ has a pseudovector coupling between the pion and the quark. The only surviving 
term for S-wave scattering is the two-pion contact interaction which comes from the 
covariant derivative on quarks; this term reproduces the Weinberg- Tomozawa result which 
was originally derived from current algebra. 

Both the scalar dielectric field and the screened gluon fields exist in the Lagrangian (^. 
The simplest treatment of them is in a mean field approximation in which we ignore the 



screened gluon fields and the pion clouds (one then needs one-gluon-exchange "corrections" 
to split the nucleon and delta masses). Pions are included perturbatively in the theory 
as we derive the bare vertex function and later dress the bare vertices by means of the 
Lippmann-Schwinger equation (|I]). 



B. Vertex Functions with the Static Baryons 

The equations of motion for the static quark and dielectric fields in MFA are: 



UYd„ -y)i = 0, (12) 

,;a„a"x + ^-^5, = o, (13) 

These are coupled, nonlinear, partial differential equations. We solve them numerically for 
x's with soliton-like behaviors. For bare N and A, all valence quarks are in IS" state, and 
so the quark wave function is assumed to be 

li^) = 4^( '^'^ ]ue-^\ (14) 



\ia-ff{r) 

where is the quark spin wave function. As usual, the scalar-isoscalar dielectric field x is 
assumed to be time-independent and spherically symmetric. 

In the same spirit as the cloudy bag model, we treat the baryons as composite, three 
quark systems, while simultaneously treating the pions as an elementary, quantized fields. 
As illustrated in Figure |l|, the interaction hadron-space Hamiltonian involving a single pion 
is then: 

< = E 4Bo. J ■ v>;m + h-c] , (15) 



a,/3 



where h.c. denotes Hermitian conjugate, a{k) is the pion annihilation operator, and the ver- 
tex function is the matrix element of the quark pseudovector current evaluated between 
transition baryon states: 
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vL(k) = / d'x 



V 



e'^-"" {Bp \qYl5rq\ B^) 



(16) 



(17) 



2/ ^{2nf2ujM 

The function u^^{k) in ([ITD is the vrNN pseudovector form factor defined in the Appendix. 
The braket in (|1^ is a spin-flavor matrix element which we reduce via the Wigner-Eckart 
theorem: 



Bp\a ■ kT\B, 



(18) 



..Sals 



IbImIbo 



where C^j^^^J^ and C^^^^^i^^ are Clebsch-Gordon coefficients and the XscMBfi are the one- 
pion coupling constants given in Table |. 

The contact interaction directly couples the initial pion and baryon to the final pion and 
baryon, Figure |l|, and takes the form: 

d^kd^k' 



HZ = Y.BlB^J (2^)3 

[{a\k') X k^aik)) ■ >V^„(k', k) + h.c. 



a, 13 
X 



(19) 



We evaluate the time and space components of the contact interaction separately: 

uj^{k) +uj^{k') 



Wjjk',k) 



A7i^P^Auj^{k)iu^{k') 
xv'l^{k',k)(Bp 



J2YUjik')YLMik) 



Br 



sf 



(20) 



W|Jk',k) 



^C!^,^^jC^^/jA'^^^v^^{k', k) ( B^ 



y: Yi,,{k')YLM{k) 



—I 



2^6 



yf 



A 



JLS 



-2V6(L||L||L)(-)-^+^+^ 



S' S 1 
L L J 



(21) 
(22) 



The spin-flavor matrix elements are 
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Br. 



Br- 



2v^ 



r,cr 



Br 



B^ 



/ 

E\SI filBhil fi^B'hi'^ 



(23) 
(24) 



where the two-pion couphng constants X^^ and are given in Table 0. In the Appendix, 
the vertex functions v'l^{k',k) and v'£^{k',k) are related to the nucleon electromagnetic 
form factors. 



C. Vertex Functions with Momentum-Projected Baryons 

The solutions of equations (p!^)-([l3|) describe quark orbits within baryons in which the 
COM may move. While this may not be too serious a concern for bound-state spectra, it 
is for scattering where the large momentum transfers make recoil effects important. We 



remove this spurious motion by a Peierls-Yoccoz [|T2| projection to form an eigenstate of 
momentum p: 

\B{p)) = ^Jdre^^-^qlqlql\C,), (25) 

where Nb is a momentum-dependent normalization constant. The operator ql creates a 
quark at point r. The ket \Cr) represents the coherent state of the dielectric field generated 
from the mean field solution 1131: 



\Cr) = exp 



|0) . (26) 



d'k^iuik)/2Mk)al,{k) 

In (P^), |0) is the vacuum state, /r(k) is the Fourier transform of the scalar dielectric field 

4b{ 



x(r), and aL(k) is the creation operator for a scalar dielectric field quanta with energy 



^{k) = yf^lb + Accordingly, the pseudovector vertex function with recoil corrections is 
evaluated between the momentum-projected baryon states in the Breit frame: 



k^ 



mri^rq{Q) 5«(|)). 



(k) = ^ ^' / ^ (27) 

After making these same projections and evaluations, the corresponding contact interaction 
Hamiltonian (ITTD takes the form: 




(i?^(-|)|g(0)7^fg(0)|i?„(f)) 



(28) 



AP^Auj^{k)u^{k') 



where q = k' — k is the momentum transfer. Detailed expressions for the vertex functions 
are given in the Appendix. 



the axial form factor of the nucleon. In particular, VN,N{k) in the A; — > limit deteremines 
qa of nucleon. We find in our calculations that qa calculated using momentum-projected 
states is larger than that calculated statically. Similar behavior is expected for other meson- 



element of the vector current, and has time and space components which are proportional 
to the charge and magnetic form factors of the baryons. 



While it is good that we have been able to solve the field equations in the mean field 
approximation, the approximations involved leave the nucleon and delta with the same 
mass. This degeneracy is removed by including one gluon exchange corrections. Since the 
scalar dielectric field x is responsible for the confinement of gluons, it is consistent to ignore 
the self- interact ion of the residual gluon. This means we drop the non-Abelian gfabcA^^^I 
term in equation (|^) for the gluon field tensor F"-''^^. Under the Coulomb gauge condition 
V • (kA") = 0, the equations for the time and space components of the gluon fields become, 



The vertex function V%a{k) is of some interest in its own right because it is related to 



baryon coupling constants. The form factor yVj^^{q) for the contact interaction is the matrix 



D. One Gluon Exchange Corrections 




(29) 




(30) 



Here j"'''^ = ^gq'~f'^X°'q is the quark current, and jj is its transverse component. We follow 



BickeboUer et al. and solve these equations for with Green's function techniques. 



As is done in bag model calculations [0, the contribution of the gluon electric energy is 



neglected. The gluon magnetic energy can then be expressed in the form: 



9 



mag 



5(0) 



a,i<j 



(31) 



where A°(i) is the gluon field due to i-th quark. The bare baryon mass is the sum of the 
expectation values of the bare Hamiltonian for a zero-momentum baryon state plus the gluon 
magnetic energy: 



mf^ = {B{0)\Hf,are\Bm+E, 



mag j 



bare 



+y{(Vx)^ + vr^} + f/(x) 



(32) 



(33) 



where vr is the conjugate momentum of the x field. When these steps are followed, the 
masses of the bare nucleon and delta can be written as 



= mgb{Ci - a^Ca), 



"^A^ = mgb{Ci + asC2), 



(34) 
(35) 



where Ci and C2 are constants determined from the mean field solutions. 



III. UNITARY 

In the present work we concentrate on ttN scattering in the energy region from threshold 
through the A resonance. Even though the vrA channel is closed for these energies, we 
include coupling to an explicit vrA channel in order to include its effect on vrN scattering 
and to permit a subsequent extension of the theory to A production. To ensure unitarity 
we iterate the tree- level diagrams using the relativistic Lippmann-Schwinger equation 
The resulting T matrices have the desired renormalized poles together with the pion-dressed 
vertex functions. In Figure |^ we show the three distinct terms which contribute in lowest 
order to the vrN and ttA potentials: 

1. Direct Born term. Figure |^. This part of the potential arises from two elementary 
vertices connected by an intermediate bare baryon state Bq. In the partial- wave basis, 
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AMP^Auj^{k)uj^{k') 

Bo E- m's'^ 

Here E is the COM energy and rn^Bo bare mass of the intermediate nucleon or delta. 

Spin and isospin conservation restricts the intermediate state to the Pn and P33 channels. 
2. Crossed Born term, Figure |^. This term arises from crossing the external pion lines 
of the pole diagram. Since the intermediate state contains two pions and one baryon, this 
modifies the propagator and so: 

r{b)LJi,,, ,^ diik'k 



VZ'\k\k) 



r]pi'u^^{k')uP^{k) 



X ^BoMB'XboM'b- — — • (37) 

Bo E - m^s'^ - uj^[k) - uj^[k') 



The spin-isospin coefficient is defined as 

\ hi Ib I W Sb I J \ 

(38) 

I Iw Ib' Ibo I I Sb' 1 Sbq I 

Recoil of the intermediate baryons are ignored in our calculations, and thus the crossed term 
contributes only to P-wave vrN scattering. 

3. Contact term, Figure |^c. The contact term arises from -^277, the two-pion vertex of 
the interaction Hamiltonian. This term gives rise to two potentials, V^J^ due to the time 
component, and V^J' due to the space component. In the partial- wave basis. 



M k)+Mk'))\^/^v'i''{k ',k) 
i7rP^4uj^{k)uj^{k') 



Vfj^'^'^k', k) = .'"^'"^'^ (39) 



V^a''''\k', k) = ^'"^'''"^'^^ ' , (40) 



where A'^^^ has been defined previously (p2l), and the spin-flavor coupling constants are 
given in Table |T|. 
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IV. RESULTS 



We have used a CDM to generate the potential terms for the relativistic Lippmann- 
Schwinger equation (|l]). The pion interaction resulting from the iteration of the potential 
by the Lippmann-Schwinger equation ensures two body unitary of T and renormalizes the 
ttBB' vertex and the bare baryon masses. The energies of the poles of T are identified 
with renormalized masses for the physical states, the residues of the poles are identified 
with renormalized coupling constants, and the momentum distributions around the poles 
are related to the form factors. 

When fitting data we always start off with a value for the glueball mass rrigb ~ 1 GeV. 
The precise value does not appear to be important, but its general magnitude does set the 
energy scale in the equations of motion. We keep the shape parameter of the self-interaction 
of dielectric field fixed at a = 24 since previous studies with the pseudoscalar version of 
the CDM [|l^ have shown that the static properties of baryons are insensitive to a for this 
size of nigb. Each time we obtain a solution of the field equations, we adjust the parameters 
nigb, as, and / to best describe the ttN phase shifts and the masses of the nucleon and 
delta. Somewhat surprisingly, we found that we can reproduce the experimental masses for 
a large range of quark masses and bag constants just by adjusting the values for rrigi, and 
as- Specifically, we find that we can vary the bag constant B^^'^ between 100 MeV and 150 
MeV, and the quark mass between 40 MeV and 120 MeV. While 40 MeV may appear small, 
the dielectric field also becomes small near the origin (the "one phase" solution in which 
the wave function has no abrupt phase transition). Consequently, the 40 MeV quark has an 
effective mass inside of a nucleon of mq/x{r ~ 0) ~ 100 MeV. 

The fitted parameters and deduced bare masses are given in Table |ITT| . The value / = 
93 MeV for the pion decay constant is the accepted value [Q, and the value / = 90 MeV is 
very close. Our value as — 0.12 for the quark-gluon coupling constant is much smaller than 



the MIT bag result, as — 2.2 or the favored value in the Friedberg-Lee model. In fact 
BickeboUer et al. |T^ have found that can vary over two orders of magnitude depending 
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on the choice of the dielectric coefficient and the parameters of the self-interaction of the 
dielectric field. In our case, a small value of compensates for our model of the dielectric 
field, K,{x) = X^) which otherwise would produce too large a color-magnetic energy. In 
addition, the pion-dressing in our model makes up for part of the and A mass difference. 
This also makes the value of a^, which measures the strength of OGE in our model, smaller 
than in those models without dressing. Our fitted value nigh — 1153MeV for the glueball 
mass is consistent with the values found in other models P|, and is in the region where 
glueball may occur jiej. Finally, we see in Table III that the effect of renormalization on the 



nucleon and delta masses is to move the bare masses down by ~ 250 MeV. Our calculation 
of scattering lengths indicates the need for some additional repulsion in our model, and we 
suspect that that repulsion would lead to less mass renormalization. 

A. Coupling Constants, Masses, and Form Factors 



As indicated in (pi]), we have deduced the CDM Hamiltonians Hi-,^ and if27r from an 
expansion of the Lagrangian in inverse powers of the pion's weak decay constant /, or 
equivalently, as a series expansion in the vrNN coupling constant f-^NN- This equivalence 
follows from the Goldberger-Trieman relation [Q: 

i- = V4^^. (41) 
2/ 

While the relation (|4l| ) relates the bare coupling constants, our study also produces renor- 
malized coupling constants. We extract the bare coupling constant /^"^tv from our final T 
by comparing the vertex function computed with the CDM to that of standard Chew-Low 
theory |[: 

v.{k)=tJ — ^g{k)T,a-k. (42) 

V 2uk 

We then extract the renormalized ttNN coupling constant Jt^nn and the renormalized ttNA 
coupling constant fnNA by making Laurent expansions of the computed T matrix around 
its poles: 
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-C/ — mjv 

fik', k;E^ m^) ^ f2^ J^-f')s^-^(^) + . . . . (44) 

This expansion permits us to deduce not only coupling constants but also the vrNN and vrNA 
form factors gnNNik) and g-rrNAik) from the momentum dependence of the T matrices in (^31) 
and p^). It is particularly interesting to compare these form factors with those calculated 
in other approaches since the form factors are related to the CDM quark wave functions 
(which should be fairly realistic), and since they are affected by the renormalization process. 

While we have already indicated that fits to the data do not determine a unique set of 
parameters, the different parameter sets are correlated. In the upper part of Figure |^ we 
show the relation between the strong coupling constant fnNN deduced from our fits and the 
bare, weak decay constant / (note the false origin). The dashed curve relates the bare values 
of fnNN as determined from the Goldberger-Trieman relation (^), or equivalently, from the 
potential directly. The solid curve relates the renormalized values determined via (^3[). We 
see that there is a similar functional dependence in each case with an ~ 15% renormalization 
effect. 

Although not shown in the figure, we have also calculated fnNA from the residue of the 
A pole and find: 

^ - 1.87. (45) 

JirNN 

This ratio is quite close to the experimental value of ~ 2, closer in fact than the SU(6) 
prediction of 6^/2/5 ~ 1.70. This clearly shows the importance of renormalization. 

In the lower part of Figure § we show the ttNN form factors deduced from the behavior 
of the T matrix near the nucleon pole. The two different curves correspond to two different 
values of the pion decay constant /. One can estimate from the falloff that the range 
R ~ 0.5fm. To compare with other models which have a monopole form factor, we assume 
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rJTlA r0.47fm / = 90MeV, 
^R = li^j9M=} (47) 
fc-o V -2k \^ 0.48 fm / = 93 MeV. 

This R is comparable with the value R ~ R\yg^g/ VW obtained from the cloudy bag model 

with a bag radius -RjQag — ^ The meson exchange models, however, use much smaller 

R ~ 0.15 fm which produces a much harder form factor.^ In practice, a monopole form 

factor falls off much slower than the form factor calculated in our model. If we replace our 

form factor with the monopole form factor defined above and use it in an actual T-matrix 

calculation, we would get poor agreement with the experimental phase shifts. Thus, in 

contrast to the meson exchange models, pion-nucleon scattering seems to require soft form 

factors like the ones obtained in our model. A similar conclusion has also been reached 

recently by others [1T7 . 



B. S and P Wave Phase Shifts 



In Figure ^ we compare the low-energy S- and P-wave phase shifts calculated with the 



CDM model to the experimental phases obtained from the said analysis |jT8[. The solid 
and dashed curves are for parameter sets with /=90 MeV and /=93 MeV respectively. In 
Table |V| a comparison of just the scattering lengths is given. We see excellent agreement 
in the P33 channel, with the energy of the delta and its width reproduced well. The small 
P-wave phases are always more of a challenge. The scattering volumes for the P31 and P13 
channels are predicted well, with the calculated energy dependence in the P31 phase very 
similar to the data. In Table |IV| we give values for the scattering lengths calculated in Born 
approximation, that is, when only the potential term in the Lippmann-Schwinger equation 
(0) is used for T. As expected, these values agree with the chiral limits: 



*We remind the reader that the monopole approximation ( |4^ ) is for comparison only and that 
this radius is determined from the behavior of the form factor in the A; ^ limit. 
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where /x is the ttN reduced mass. Comparison of the Born results to experiment and to our 
full calculation indicates that renormalization is large, and that some additional repulsion 
is needed. 

The Pii channel is a particular problem for simple models like ours. Increasing the agree- 
ment with the P33 phases decreases the agreement with the Pn phases. We have assumed 
that the Roper(1440) is a (ls)^(2s) three-quark state, and included it as an intermediate 
state for both s- and u-channel diagrams. The results show that the elementary Roper does 
provide the attraction needed for the phase shift to change sign from negative to positive, 
but does not do much to improve the fit at lower energies. In addition, we find that the Pn 
amplitude does not have the desired pole at the physical Roper resonance unless we intro- 



duce extra, free parameters. Indeed, the calculations of Elsey and Afnan and Pearce 
and Afnan pO|] indicate that the coupling of three-body (vrvrN) channels is important, as 
are some 1//^ diagrams not present in our calculation. These should be included in the 
continuation of the present work. 

As is true for the small P-wave phases, the predicted S-wave phase shift shown in Figure |^ 
show the need for additional repulsion. Cooper and Jennings |^T]] indicate that agreement 
may require the inclusion of certain 1//^ graphs which cancel at m.^ = but produce 
repulsion for m^r > 0. 



V. SUMMARY AND CONCLUSION 

We have studied low-energy vrN elastic scattering in a chiral version of the color dielectric 
model. Our calculations are similar to those done with the cloudy bag model 0. Our major 
improvement is a more natural and realistic inclusion of quark confinement, the inclusion 
of baryon-recoil effects, and the coupling of the vrN and vrA channels. The CDM is used 
to derive effective potentials containing terms up to second order in the vrNN coupling, 
and these potentials are used as input to relativistic, coupled-channel Lippmann-Schwinger 
equations. The model reproduces the physical nucleon and delta masses, and then predicts 
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the ttNN and vrNA coupling constants and form factors, and low-energy vrN scattering. 

We find much success for such a simple model. Specifically, we can simultaneously obtain 
the correct renormalized masses for the nucleon and delta, the correct delta width, the vrNN 
coupling constant to within 5% of it experimental value, the ratio of the renormalized vrNA 
to vrNN coupling constants which is better than the SU(6) value, excellent agreement with 
the P33 phase shifts from threshold through the delta resonance energy, good agreement with 
the Pi3 and P31 scattering volumes, and good agreement with the energy dependence of the 
P31 phase shift for the first 300 MeV of kinetic energy. Better agreement with the S-wave 
and small P-wave phases requires additional repulsion in the model, possibly obtained by 
one-loop corrections. The model could also be improved by obtaining solutions to the field 
equations beyond those of the mean field approximation and by the inclusion of additional 
dielectric fields p2[ . 



APPENDIX: FORM FACTORS FOR ttN SCATTERING 

We give here the complete form factors for the (7rN,7rA) system. Since the A differs from 
the only in spin-isospin orientations, they have the same space wave function and thus 
the same form factors. The pseudovector vertex functions with the momentum projected 
states evaluated in the Breit frame are: 



vL(k) = _ 



2fJ2u;„{k) 



2fpu^{k)NB,NB^ 
u^^{k)(Bp a -kr . (Al) 



Where Nb^ = y \Ba{2)\Ba{2)/ is a momentum-dependent normalization constant for a pro- 
jected baryon state, and |Pq)x = it.li.Q^ |C'x) is a localized baryon state. The pseudovector 
form factor has the form, u^^{k) = N{k)/D{k), with 
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N{k) = J z^dzN^{z)N^{z) J 



X i 9+9-jo{kr) 



/-/+ 
3r_r+ 



bo(^?^) + 4j2(/i;r)]r^ - \jo{kr) + Aj2{kr)P2{f ■ z)] — 



D{k) = j z'dzjo{kz/2)N^{z)N^{z), 
N,(z) = / d\ + ^ (r^ - J 

N^{z) =exp 27r k'^dkuj{k) \fo{k)f jo{kz) 



(A2) 
(A3) 
(A4) 
(A5) 



In the above equations, we have defined r± = |r ± z/2|, g± — g{r±), and f± — f{r±). The 
function fo{k) is the Fourier transform of the mean field solution of the classical dielectric 
field. Without the COM correction, g± — g{r), f± = f{r), and the ttNN pseudovector form 
factor is reduced to the static limit: 



„^^(A;) = N'i / r'dr 



9{rY - 



2 f\r) 



o — J2{kr) 



The vertex functions for the contact interaction are 

1 



>vL(q) 



SP^uj,{k)uj,{k') 

1 



X I d^dye-"^<'^+yy\{Bp\mrMQ)\Ba)^ 
1 

8Py/u;^{k)u;^{k>) 

u^^{q){Bp\f\B^f , time, 
u'^^iq) {Bp\ -ax qf|Sc,)^-^ , space. 
The form factors for the time and space components are respectively: 



X 



u^'{q) 



N^^iq) 



D{q) ' D{q) ' 

N^^{q) = j z^dzNl{z)N^{z) j d^ rjo{qr) 



X 



'g^g +(r'--]^ 



(A6) 



(A7) 



(A8) 



(A9) 
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2 ( 9+f- , ^ 



X 



+ 



qr 

r-zf g+f. 



r+ 



(AlO) 



Their corresponding static expressions are: 



u 



CT 



= Jdrr'jo{qr)[g'{r) + f{r) 



drr 



g{r)f{r) 



(All) 
(A12) 



Finally, the relationships between the electromagnetic form factors and the partial-wave 
decomposed radial potentials are: 



v?Hk', k)^^ dx [1 - P^ix)] ucsiq). 



(A13) 
(A14) 
(A15) 
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TABLES 








TABLE L One-pion spin-flavor coupling constants. 










ttN 




ttA 


N 
A 


5 




4^2 
5 


TABLE 11. Two-pioii si)iii-(lav()r c()ui")liiig cousiMiiis. 






1=1 /=3 /=1 
^2 ^2 ^2 

ttN ttA ttN ttA ttN ttA 


ttN 


2 

ttA 


ttN 
ttA 


-1 i -f 

5 .1 .272 

2 J D 


5 
6 

275 
3 


2a/5 
3 

a/Io 

3 


TABLE in. Fitted parameters and deduced bare masses for the CDM. 


/ (MeV) 


as (MeV) mj^^ (MeV) 




m^^ (MeV) 


90 
93 


0.112 1178 1239 
0.126 1127 1171 




1489 
1439 
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TABLE IV. S and P wave scattering lengths in fm where (kaa) '^^ = tanSa- 



L2I2J 


acDM 
(/=90 MeV) 


CtCDM 

(/=93 MeV) 


(^Born 


aCBM 


Exp. 


Su 


+0.398 


+0.342 


0.237 


+0.42 


+0.243 




-0.098 


-0.093 


-0.119 


-0.07 


-0.130 


Pll 


-0.892 


-0.860 


0.954 




-0.569 


P3I 


-0.412 


-0.420 


-0.426 




-0.483 


Pl3 


-0.355 


-0.366 


-0.357 




-0.375 




+0.713 


+0.712 


0.583 




+0.839 
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FIGURES 



FIG. 1. The linear and quadratic term in 1// contributing to the interaction Hamiltonian. 

FIG. 2. The three lowest-order terms of the potential contributing to ttN scattering. 

FIG. 3. {Top) The relation between the ttNN coupling constant and the pion decay constant 
/. {Bottom) The ttNN form factor for two different values of the pion decay constant /. 

FIG. 4. The vrNN form factors deduced from the behavior of the T matrix near the nucleon pole. 

FIG. 5. The S- and P-wave vrN phase shifts. The solid and dashed curves are for parameter 
sets with /=90 and /=93 respectively, and the dots are the experimental phase shifts. 
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